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MOLECULAR ATTENUATION AND PHASE DISPERSION 
BETWEEN 40 AND 140 GHz 
FOR PATH MODELS FROM DIFFERENT ALTITUDES 

Hans J. Liiebe and W. M. Welch 


ABSTRACT 


Radio wave propagation in the 40 to 140 GHz band 
through the first hundred kilometers of the atmosphere 
is strongly influenced by the microwave spectrum of 
oxygen (O 2 -MS). A unified treatment of molecular 
attenuation and phase dispersion is formulated. Results 
of molecular physics are translated into frequency, 
temperature, pressure, and magnetic field dependencies 
of a complex refractive index. The intensity distribution 
of the O 2 -MS undergoes several changes with increasing 
altitude: when h < 10 km, all lines, but one at 119 GHz, 
are merged to a continuum spectrum under the influence 
of pre s sure -broadening; when h > 30 km, a line spectrum 
with isolated Lrorentzians is displayed; when h > 40 km, 
Zeeman-splitting of each line occurs due to the influence 
of the earth^ s magnetic field; for h > 60 km, a Voigt pro- 
file governs the transition to a Gaussian line shape and 
eventually the Doppler -broadened line spectrum vanishes. 
The influence of water vapor is discussed separately. 

Attenuation and dispersion rates for path models are 
evaluated by computer routines. Examples of computer 
plots are given as a function of altitude for homogeneous, 
zenith, and tangential path geometries. Molecular reso- 
nances of minor atmospheric gases are discussed briefly, 
as in the noise which originates from the O^-MS. 

Keywords: Atmospheric transmissivity; atmospheric 
millimeter wave propagation; complex refractive 
index of air; oxygen microwave spectrum. 


^The authors are with the Institute for Telecommunication Sciences, 
Office of Telecommunications, U. S. Department of Commerce, 
Boulder, Colorado 80302. 
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1. INTRODUCTION 

Much of the future growth of telecommunications will occur by 
extending the radio spectrum into the millimeter wavelength range (30 to 
300 GHz). As shown in a recent review by Thompson et al. (1972), the 
successful exploitation of this largely unused portion of the spectrum 
hinges on the accurate and complete account of all propagation effects. 
The greatest obstacles to the use of such high frequencies are attenuation 
and phase dispersion caused permanently by the clear air and inter- 
mittently by rain and clouds. This report addresses the use of the fre- 
quency spectrum 40 to 140 GHz for radio propagation through the clear 
atmosphere and discusses in detail ultimate limitations due to air mole- 
cules of which oxygen is a major species. The microwave spectrum of 
oxygen, abbreviated O^-MS, dominates the transfer properties of air 
throughout the specified frequency range, while molecular effects of 

water vapor influence them to a lesser extent. Roughly forty lines of the 
O^-MS contribute to the atmospheric spectrum. Each line has individual 
pressure and temperature dependencie s, and the overlapping of lines 
below 30 kilometers changes the structure of the spectrum drastically 
with altitude. 

Of prime concern to system engineers are signal degradations 
(amplitude and phase distortions; interferences) and the minimum 
transmitted pulse length supported by the propagation medium. These 
concerns lead in Section 2 to the formulation of a complex transfer 
function T(a, 0) for a well -mixed, homogeneous atmosphere. Spectros- 
copic knowledge of the O^-MS is translated into engineering terms of 
attenuation a , phase -delay 0^, and dispersion A0 rates for atmospheric 
conditions up to altitudes of 100 kilometers. In Section 3 these results 
are applied to calculate transfer properties of path models through the 
U. S. Standard Atmosphere (1962). Cumulative attenuation A and phase 
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dispersion AT along slanted ray paths through the total (inhomogeneous) 
atmosphere are evaluated by computer calculations. An example is 
shown in Figure 1 for one-way zenith path attenuation between 49 and 
72 GHz from different initial altitudes. A graphical presentation is the 
best introduction to atmospheric transmission and shielding properties 
due to the O^-MS. 

As part of spectrum management, a first step was taken to 
allocate bands for specific applications through an international agree- 
ment (ITU, 1971), For the first time several bands have been reserved 
exclusively for passive systems, A clear understanding of the O^-MS 
is required to promote the most efficient use of the millimeter wave 
spectrum which, just as its lower frequency counterpart, is a saturable 
resource. Reliable knowledge of the transfer function t provides the 
basis for designing systems which take advantage of unique molecular 
transmission and emission properties associated with the relatively 
stable dry part of the atmosphere. The O^-MS affords interference 
resistance and transmission security for broadband radio communica- 
tions. Optimum system performances can be assured by tradeoff 
studies between environmental degradations and component limitations. 

The following key phrases give examples of telecommunication 
applications which incorporate O^-MS properties: 

(a) Shielding of broadband satellite -to -satellite links from 
ground interferences (e.g., Crane, 1971). 

(b) Secure high -altitude communication in the valleys of the 
O^-MS with special channel characteristics (balance 
between transmission and shielding, phase -dispersive 
pass-band, phase keying, etc.)- 

(c) Restricted range, variable range, and local distribution 
broadband communication systems (e.g., Murray, 1971). 
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(d) Aircraft collision avoidance by means of «==«60 GHz beacons 
and "threshold'* receivers (M. L». Meeks, private communi- 
cation, 1969)* 

(e) Restrictions of radar range (Blake, 1972). 

(f) Occultation experiments between orbiting satellites for 
studies of atmospheric structures. 

(g) Absolute differential phase measurements between frequency 
pairs on either side of the O^-MS over line -of-sight paths 
for correction of radio distance (MITRE, 1965; Thompson, 
1968; see Section 3. 1. 1). 

More applications will almost certainly follow if the basics are better 
known. 

Active systems may eliminate desirable molecular properties 

(e. g. , shielding) when operating at extremely high power densities. A 

very rough estimate indicates that, for the onset of saturation at an 

2 2 

isolated line center, the power density has to exceed 10*p [W/(cmtorr) ]. 
The nonlinear power saturation effects, however, are not treated any 
further in this report. 

The following applications are based on frequency -selective 
emission originating from the O^-MS: 

(a) Remote sensing of atmospheric temperature structures 
(Meeks and Lilley, 1963; Lenoir, 1968; Wilheit, 1969; 
Waters, 1970; Westwater, 1970). Ground-based strato- 
spheric temperature soundings (Waters, 1973) and the 
NIMBUS-5 (launched Dec. 1972) satellite (Sabatini, 1972), 
both operate radiometers within the O -MS band. The con- 
version from measured radio noise to meteorological quan- 
tities requires accurate knowledge of oxygen absorption. 

This application in particular has fostered theoretical and 
experimental work on the O^-MS (see Sect. 3. 2. 4). 
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(b) Stratospheric navigation aids such as velocity and altitude 

/ 

indicators (Mardon, 1969). 

(c) Horizon sensor for satellite navigation (Guidice, 1971; 
Grauling, 1972), 

(d) Radiometric detection of strong temperature gradients in 
the atmosphere for the identification of clear air turbulences 
(Haroules and Brown, 1969). 

(e) Earth magnetic field strength measurements from satellites 
by means of the normal Zeeman effect {Sect, 2, 3, 3 and 
3,1, 1), 

(f) Line emission (e. g. 119 GHz) used as a tracer in 
biological re search, for example, to study oxygen exchange 
rates. 

Radio astronomy efforts which are searching for millimeter wavelength 
emissions from exterrestrial sources are interested in accurately 
defining the boundaries of the atmospheric O^-MS, which masks their 
useful frequency bands. 

This report shall provide baseline data on unique atmospheric 
transfer properties for a largely unused portion of the spectrum. 
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Zenith Attenuation 



Figure 1. Millimeter-wave attenuation in the 49 to 72 GHz band 
due to atmospheric oxygen resonances for zenith paths 
through the U. S. Standard Atmosphere 1962 from dif- 
ferent initial heights h^ to outer space (see also Section 
3.2). Allocated band utilizations are: A(1 -5) -active 
telecommunications; P(l-3) - passive systems only 
(I. T. U. , 1971; Klass, 1971). 
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Z. THEORY 

2. 1 Atmospheric Transfer Function 

Signal analysis of millimeter waves propagating through a well 
mixed atmosphere is facilitated by treating the medium as a linear 
system and applying the transfer function concept, (Morgan and Ekdahl, 
1966). The central topic is the development of a microwave transfer 
function for moist air as reliably as can be at the present state of know- 
ledge. To accomplish this, one has to track down macroscopic pro- 
perties of the propagation medium to their molecular origin and be con- 
fronted with a combination of classical and quantum mechanical theories 
which describe the interaction between radiation and gas molecules. A 
rigorous test of these theories is still wanting since it can not be pro- 
vided by the small amount of inaccurate and sometimes conflicting 
experimental data. The complex transfer function is defined by 


1 ^ ~ ^ -( 20^og e U 2nv La/ c] (1) 

where is the initial amplitude and E is the amplitude and phase 
received after the geometric path length L for free space propagation; 
j =vtt: Power attenuation rate a and phase -delay rate 0 are the well 
known propagation parameters 


a = ( 10 log e) (4tt v/ c)n^^ [dB/ km] , 
0 = (2 tt v/ c)n^ [rad/ km] 


( 2 ) 


where V is the microwave frequency, c is the speed of light, and 
10 log e 4. 34294 . The dimensionless macroscopic measure of the 
interaction between radiation and gas molecules is a complex refractive 
index 
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n = n' - jn” , (2a) 

where 

n' = n +^ri(v) is the Refraction Spectrum and 
o 

n” (v) is the Extinction Spectrum. 


Both spectra are interrelated by Kramers-Kronig integral equations, 
which are obeyed by the shape functions chosen (Eqs. 17, 23, 28, 
Table 8). This report is concerned with theoretical and analytical 
evaluations of the refractive index n for air over the frequency band, 

V « 40 to 140 GHz. The environmental conditions to be considered for 
altitudes h = 0 to 100 km are: 


Dry air pressure 

Mixing ratio of dry air 
components 

Temperature 

Water Vapor pressure 

Earth's magnetic field 
strength 


P(h) 


[ see below] 


r, =const. 
k 

T(h) 

P«<T. h) 

H(h and geo- 
magnetic 
coord. ) 


[see Table l] 

[330 to 180“K] 

[O to 50 torr]* 

[O. 2 to 0. 7 gauss]* 


(3) 


The field strength H decreases between h = 0 and 100 km by 4. 8 percent. 


Pressure changes versus altitude as shown (U. S. Std. Atm. 62): 


h [km] 0 100 20 30 40 

p [torr]* 760 199 41.5 9-0 2.2 


h[km] 50 60 70 80 100 

pCmtorr] 600 170 41 7.8 0.23 


^The units [torr] and [gauss] do not comply with recommended S. I. 
units. The conversions are 1 torr = 133. 322 pascal (N/ m^) and 
1 gauss = 10-^ Wb/m3. 
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The solution of equations (1) and (2) leads to the extinction and 
refraction spectra of oxygen and water vapor being defined by 


and by 


n' 


= n 


S S. F! 
. 1 1 
1 


+ An 


n 


An 


w 


S S. FV + n” + n" 

. 1 1 n w 

1 


(4) 


where 

n^ is the frequency -independent refraction of air (see Sect. 2. 2); 
i is the label for lines of the Microwave Spectrum 
(see Sect. 2. 3. 1); 

Sj^ are the individual line strengths (see Sect. 2. 3. 2 and 3); 

F^, F." are the shape functions (real and imaginary parts) for 
individual lines (see Sect. 2.4); 
n” and An indicates the nonre sonant oxygen spectrum (see 
Sect. 2. 4. 2); and 

n” and An indicates contributions due to water vapor (see 
w w ^ 

Sect. 2, 5). 

2, 2 Frequency-Independent Phase Delay 

The non -dispersive phase delay of a millimeter wave propagating 
through the gaseous atmosphere is commonly expressed in N -units, 

N = (n -1) 10^ = [pR; + p • R* (T)]/T. [ppm] (5) 

The corresponding phase delay in engineering terms is 

0 =(2rrv/c)-n , [rad/ km] 

o o 

^Pressure proportionality is indicated in general for all appropriate 
parameters by a superscript ^ . 
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Table 1. Relative Concentration r of Dry Air Components (near sea 
level) and Refractivity R Measured for Different Air 
Components at v = 61 156 MHz, p 1 20 torr, T = 280 ° K. 
(Liebe and Welch, 1972). 


k 

Component 


r 

R°(280) 

R°/T 




ppm by vol 

ppm/torr 

ppm^^K/torr 

1 

Oxygen 

03 

209 460^®^ 

0. 3427 

95 . 95 

2 

Nitrogen 

Ns 

780 840 

0. 3771 

105. 60 

3 

Argon 

Ar 

9 340 

0. 3562* 

99. 73* 

4 

Carbon Dioxide 

CO 8 

314 

0. 6356 

178. 0 

5 

Neon 

Ne 

18.2 

0. 08659 

24. 25 

6 

Helium 

He 

5.2 

0. 0450 

12.6 

7 

Methane 

CH 4 

2 

0. 5672 

158. 8 

8 

Krypton 

Kr 

1. 1 

0. 5428 

152. 0 

9 

Nitrous Oxide 

NgO 

0.5 

0. 7138 

199. 9 


Dry Air 


10® 

0. 3696 

103.54 


Water Vapor 

Hg 0 

variable 

6. 712 

1 879(280° K) 


(a) Machta and Hughes (1970). 

* Reference value; all measured R° are relative to the R° (Ar) - 
value given by Newell and Baird (1965). 
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The refractivity values for dry air and water vapor are 
R * = 2 r R “ = 103. 5 [ppm*K/ torr] 

Q , K. K 

and ■' 

R* = 95. 5 + 499500/ T [ppm*K/torr] , 
w 

which were measured 'for the air components at 6l. 156 GHz as listed in 
Table 1^ Similar results have been reported by Boudouris (1963) and 
Newell and Baird (1965), 

2. 3 The Microwave Spectrum of Oxygen 

The following simplified explanation can be given as to the nature 
of the O^-MS. The magnetic dipole moment of the molecule couples 
to the microwave radiation field which excites quantized changes 
(molecular resonances) between fine structure levels of the rotational 
energy states. At these resonances the molecules momentarily store 
(phase dispersion), absorb (attenuation), and randomly re -radiate 
(incoherent noise) microwave energy. 

Solutions to problems of atmospheric millimeter wave trans- 
mission (Eq. 4) require detailed knowledge of the Microwave Spec- 
trum. Parameters of the O^-MS which originate from properties of 
^isolated” molecules are the line center frequencies and the strength 
values S. 

2.3. 1 Line Center Freque^ncies of the O^-Microwave Spectrum 

The O fine structure microwave spectrum obeys the selection 
rule: AN=0, AJ = ± 1 (N, J are the rotational and total angular momentum 
quantum numbers), and the transition frequencies can be calculated 
by means of (Welch and Mizushima, 1972) 
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(6) 


(N-) = \ + ytJL - + (X + I |i ) (a® + a + 2) 

O 0^0 Xl^l 

± {f (a) - (2 a + 1) [ P (a) + J X^ - n^] } , 

where a = N±l, N = 1, 3, . . . (only odd values are allowed), 
f (a) = V + a(a+ 1)I>2 , 

C = (2a+ 1) [P(a) (a®+ a + 4)] 

- [Xo+ j X^ (7a2+ 7a+ 4)] / (2a+ 1) 

D = 2[X^+ X^ (a^+ a+ 1)] /(2a+ 1) 

p (a) = B + 2B (a^ + a + 1) + B (3a^+ 6 a® + 13 a® + 10 a + 4) . 

° 1 2 

The line designation N denotes the transitions J = N ± 1 — J' = N. 
The most recent values of the molecular constants are given by Welch 
and Mizushima (1972) and they are in megahertz: 

Rotational: B^ = 43100. 518 (3), Bj = -0.14492 (9), 

Ba = 1. 57 (11) • 10"’ 

Spin-spin coupling: X^ = 59501. 342 (7), 

Xi = 5. 847 (3) . 10"® . 

Spin -rotational coupling: \ jl ^ = -Z52. 5865 (10), 

Hi = -2.464 (20) . lo"^ . 
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For all atmospheric transmission problems, it is sufficient to 

consider i - 44 individual lines. The known values of are given in 

o 

Table Z. Most of the center frequencies have been measured with high 
accuracy, and the few missing ones are predicted by equation (6). 

These predictions, for example, gmded atmospheric emission measure- 
ments by Waters (1973) whereby the agreement was excellent (^ee 
Table Z). 

2. 3, Z Line Strengths of the O -Microwave Spectrum 

Ct 

For practical applications it is helpful to reduce molecular line 
shape theory and define a line strength parameter S (Liebe, 1969a) as 
a multiplier to a Lorentzian shape function (Sect. Z. 4. 1). Pressure and 
temperature dependences of the strength for an individual line of the 
O^-MS are given by 

S= S.o (T) §p . [Hz] (7) 

. * 

where =0. 20946 (Table 1), and § = 0. 9976, portion of molecules of 
mass . 

The strength parameter S° was derived from the equation des- 
cribing the transition probabilities. Assuming for the magnetic field 
strength H = 0 (H ^ 0, see Section 2, 3- 3) leads to the expression 
(VeinVleck, 1947; Morgan and Ekdahl, 1966; LeFande, 1968; Wilheit, 
1969) 

S* = . K . f(N) . f(T) , [Hz/torr] (8) 

where K = 69. 50125. 10“^.k(J) [ torr] and k(J) is given in Table 3 
and discussed below. 


The isotopic abundance of percent 

re spectively. 
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The intensity functions (Ivlatrix elements) are 


= (2N®+ 3N) / (N+ 1) 
f(N") = (2N®+ N- 1) /N . 

3 

The temperature function is in units of [1/ “K ] 

f(T) = -[T“^exp[-2. 06858. N(N+1)/Tj-[1- (hv^ /2 k T)] 

The error of the approximation [l - / ZkTj^ 1 is less than 1 per- 

cent for all values of interest (see Fig. 3)* 

The approximation made in calculating the matrix element of the 
magnetic dipole moment operator was that is a pure Hund^ s case 
molecule. The fact that this is not exactly true was accounted for 
by a correction factor k(J) ^ 1, which is a mixing coefficient quantifying 
the purity of a particular state (Zimmerer and Mizushima, 1961; Wilheit, 
1969). The k(J) values given in Table 3 were obtained by diagonalizing 

the Hamiltonian matrix with the appropriate wave functions. Equation 

i 

(8) was applied to compute at T = 300 up to N =43, and the results 
are listed in Table 3. Figure Z displays the strength value distribution 
versus the frequency positions for ZOO® and 300 ®K. 

The temperature dependence of the strength parameter can be 
expressed relative to the reference temperature, T^ = 300K, by 

S“ (T) = S.» (300) . ^(N) , (11) 

with 

^(N) = (300/ T)^. exp j^-6, 89526.10^^- N(N+ 1)[(300/ T)-l]} . (12) 

The temperature function can be approximated for a small temperature 
interval (A T« 40 °K) by 

^«(T^/T)^^^^ . (13) 


( 9 ) 


( 10 ) 
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Each N-pair of lines possesses its own temperature dependence <>3^(N). 
The values of are given in graphical form in Figure Z for atmospheric 
temperatures, and numerical values are listed in Table 4 which also 
shows temperature exponents w(N) of equation (13) for := 280 ± 20*K. 
The temperature function '0' indicates the population of energy states 
relative to 300 ®K. At low temperatures many of the lines with high 
N-numbers may be neglected. 


2, 3. 3 Zeeman-Splitting of the Lanes in the Microwave Spectrum 

The O^-MS is greatly complicated by Zeeman -splitting of each 
fine -structure line into many components under the influence of a per- 
manent magnetic field such as the earth' s. The Zeeman effect, which 
originates from the interaction of the molecular magnetic moment with 
a static magnetic field of strength H, removes the spatial M-degeneracy 
of the J levels, splitting each level into (2J + 1) sub-levels (M-magnetic 
quantum nximber). The complete Zeeman Hamiltonian for was given 
by Tinkham and Strandberg (1954). In the weak field case (H < 1 gauss), 
however, it is sufficient to retain only the term involving the interaction 
of the molecular electronic spin moment with the field. Matrix elements 
of the spin moment, evaluated in the Hund' s (case "b") representation, 
are given by Evenson et al.(1968) and may be used in a perturbation 
calculation to obtain the Zeeman transition frequencies. 

The selection rule on AM (change in magnetic quantum number) 
for Zeeman transitions is derived from assumptions about the orientation 
between H and the microwave magnetic field component perpendicular 
to the direction of propagation. By a judicious choice of the orientation 
between H and , one can selectively excite two different types of 
Zeeman transitions: 
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(1) Linearly polarized radiation ^^induces^’ 


AM=0, called TT transitions when H and are perpendicular, 
AM = ± 1, called a transitions when H and ^ are parallel, and 
in both cases H is perpendicular to the propagation direction. 

(2) Circular polarized radiation ^’induces" 

AM= + 1 or transitions for right-handed polarization, 

AM = - 1 or a transitions for left-handed polarization, and 
in both cases H is parallel to the direction of propagation. 

The simplest pattern, the normal Zeeman effect, is exhibited by 
the 1^ and 1 lines (see Table 5), each splitting into three components 
whereby only the two a components shift away from the center (Hill and 
Gordy, 1953), The anomalous Zeeman effect splits the lines N = 3 , 

± , ❖ X 

5 , , , * into 3(2 J + 1) components where J is the smaller of the two 
J values involved in the transition. 

The frequency of each Zeeman component is determined by 


Vo = Vq + T](M, N) . ZH [MHz] . 


The shift is proportional to the external field strength H, and r[ (Table 5) 
is less than one. The factor of proportionality for ^^Icrowave lines 

is Z = 2.8026 [MHz/ gauss]. The earth’ s magnetic field strength, 

H = 0. 2 to 0. 7 gauss, produces frequency shifts for the Zeeman com- 
ponents of up to 2 MHz. 

The line strength of each Zeeman component is given by 

S° = . K .f (M, N) . f(T) , [Hz/torr] (15) 

z o 

analogous to equation (8). The intensity function f(M, N) is listed for 
the different Zeeman components in Table 5. Summing over polarizations 
and over all M - degeneracies of the final state J’ yields f(N). 

Examples for atmospheric conditions are discussed in Section 3. 1. 1. 
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Table 2, Observed and Calculated Tine Center Frequencies v ^ 
of the ^^0 2 Microwave Spectrum, in MHz* 


L 

Line 

Observed 

Ref. 

Calculated 
(Eq. 6)'") 

i 

Line 

N' 

Observed 

Ref. 

Calculated 
(Eq. 6)*> 

1 

1‘^ 

56264. 766(20) 

a 

56264. 758 

2 

1 - 

118750. 343(10) 

b 

1 J8750. 330 



56264. 778(10) 

b 








3+ 

58446. 580(10) 

c 

58446. 580 


3- 

62486. 255(10) 

c 

62486. 267 



58446. 600(10) 

b 




62486. 225(10) 

b 


5 

5 + 

59590. 978(10) 

c 

59590. 979 


5" 

60306. 044(1 0) 

c 

60306. 055 


7 + 

60434. 776(10) 

c 

60434, 778 


7- 

591 64. 215(10) 

c 

59164. 21 1 



60434. 788(10) 

f 








9 + 

61 150. 570(1 0) 

c 

61150. 567 

10 

9“ 

5832 3. 885(10) 

c 

5832 3. 883 



61 1 50. 565( 5) 

f 








11 + 

61800. 169(20) 

a 

61800. 167 


11- 

57611. 4 ( 2) 

d 

57612. 49Z 



61800. 155(10) 

c 








1 3+ 

62 411. 22 3(10) 

c 

62411 . 234 


13- 

56968. 180(20) 

a 

56968. 214 

16 

15 + 

62996. 6 { 2) 

d 

62997. 999 


15- 

5636 3. 393(20) 

a 

56363. 397 


17 + 

63568. 520(10) 

c 

63568. 542 


17- 

55783. 819(20) 

a 

55783. 805 


19 + 

64127. 777(20) 

a 

64127. 790 

20 

19“ 

55221. 372(20) 

a 

55221. 362 


21 + 

64678. 2 ( 2) 

d 

64678. 920 


21 ' 

54671. 145(20) 

a 

5467L 141 


2 3+ 

65224. 12 (20) 

c 

65224. 076 


23- 

54129. 4 ( 4) 

d 

54129. 926 








541 30.2 ( 5) 

e 


25 

25 + 

65764. 744(20) 

a 

65764. 760 


25- 

63599. 4 ( 8) 

d 

53595. 682 








53595. 9 ( 2) 

e 



27 + 



66302. 06 

1 

27- 

53066. 9 ( 2) 

e 

5 3066. 802 


29 + 



66836. 77 

30 1 

29' 

52542. 4 ( 2) 

e 

52542. 23 


31 + 



67369.51 


31- 

52021. 4 ( 5) 

e 

52021. 17 


33+ 



67900.73 


33- 



51503. 02 

35 

35 + 



68430. 8 


36- 



50987. 3 


37+ 



68960. 1 


37- 



50473. 6 


39 + 



69488. 7 

40 

39- 



49961. 8 


41 + 



70016. 9 


41- 



49451. 4 

43 

43+ 



70544. 9 

44 

43- 



48942. 4 


( ) Uncertainty of the last digits given 


*^^See also Welch, and Mizushima (1972) 


References - a. West and Mizushima (1966) 

b. McKnight and Gordy (1968) 

c. Zimmerer and Mizushima (1961) 

d. Mizushima and Hill (1954) 

e. Waters (1973) 

f . Liiebe and Welch (1972) 
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Table 3, Lrine Strength S® (300) at T = 300^K (Eq» 8) and Strti^ngth Correction k(J) 
for the Microwave Spectrum, 


i 

Line 
N + 

k(J) 

(300) 

[ Hz/torr] 

i 

Line 

N“ 

k(J) 

S° (300) 
[Hz/ torr] 

1 

r 

0. 9808 

0. 348 69 

r 

2 

1" 

1 

0. 597 25 


3-^ 

0. 9941 

0.925 05 


3- 

0. 9808 

0. 963 36 

5 

5^ 

0.9972 

1 . 341 02 


5- 

0.9941 

1. 348 70 


7^ 

0. 9984 

1. 562 63 


7- 

0.9972 

1. 526 31 


9'" 

0. 9989 

1. 589 94 

10 

9- 

0. 9984 

1.515 04 


11 + 

0. 9992 

1. 458 82 


11- 

0. 9989 

n 359 54 


13 + 

0. 9997 

1.227 20 


13- 

0.9992 

1. 120 44 

15 

15 + 

1 

0. 954 01 


15- 

0. 9997 

0. 853 87 


17 + 

1 

0. 689 76 


17- 

1 

0. 605 61 


19 + 

1 

0. 465 60 

20 

19- 

1 

0. 401 19 


21 + 

1 

0. 294 21 


21- 

1 

0. 248 87 


23 + 

1 

0. 174 37 


23- 

1 

0. 144 84 

25 

25 + 

1 

0. 097 074 


25- 

1 

0. 079 186 


27 + 

1 

0. 05 0 82 0 


27- 

1 

0. 040 717 


29 + 

1 

0.025 041 

30 

29" 

1 

0. 019 707 


31 + 

1 

0. on 621 


31- 

1 

0, 008 985 


33+ 

1 

0. 005 083 


33- 

1 

0. 003 860 

35 

35 + 

1 

0. 002 096 


35- 

1 

0. 001 564 


37 + 

1 

0. 000 815 


37- 

1 

0. 000 598 

39 

39 + 

1 

0. 000 299 

40 1 

39- 

1 

0. 000 216 


41 + 

1 

0. 000104 


41“ 

1 

0. 000 073 

43 

43 + 

1 

0. 000 034 

44 

43 - 

1 

0, 000 024 


Line Strength S° , Hz/torr 



Figure 2. Line strength parameter for the O 2 -MS at 200°K and 300'’K. 
Numerical values at 300 °K are listed in Table 3. 
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Table 4, Temperature Dependence V of the Line Strength S®(300) of the O Microwave Spectrum 


LINE 

N± 

180"K 

200 "K 

220“K 

Tempera 
240 "K 

ture T 
260 'K 

280 °K 

300°K 

320 '’K 

340"K 

w 

(280i 

20»K) 

1 

4. 5873 

3. 3518 

2. 5230 

1. 9464 

1, 5329 

1, 2287 

1 

0. 8247 

0. 6881 

2. 99 

3 

4. 3812 

3. 2382 

2,4605 

1. 9131 

1. 5168 

1, 2227 

1 

0. 8283 

0. 6937 

2, 91 

5 

4. 0333 

3. 0434 

2. 3519 

1. 8547 

1. 4881 

1. 2119 

1 

0. 8347 

0. 7039 

2. 78 

7 

3, 5789 

2. 7824 

2. 2035 

1. 7734 

1. 4476 

1. 1965 

1 

0. 8441 

0. 7189 

2.59 

9 

3. 0611 

2. 4747 

2. 0234 

1. 6724 

1. 3963 

1, 1766 

1 

0. 8566 

0. 7390 

2. 33 

1 1 

2. 5236 

2. 1411 

1. 8212 

1. 5556 

1. 3355 

1. 1525 

1 

0. 8722 

0. 7646 

2. 02 

13 

2. 0054 

1. 8021 

1 . 6066 

1.4272 

1 . 2665 

1. 1245 

1 

0. 8912 

0. 7962 

1 . 65 

15 

1. 5361 

1. 4755 

1. 3892 

1. 2914 

1. 1909 

1. 0928 

1 

0. 9138 

0. 8346 

1.22 

17 

1. 1341 

1. 1752 

1. 1773 

1, 1525 

1. 1104 

1. 0579 

1 

0. 9401 

0. 8805 

. 73 

19 

0. 8071 

0. 9106 

0. 9779 

1. 0145 

1. 0265 

1. 0200 

1 

0. 9706 

0. 9350 

+ . 18 

21 

0. 5536 

0. 6863 

0. 7962 

0. 8808 

0. 9410 

0. 9796 ' 

1 

1. 0055 

0, 9993 

- . 42 

23 

0. 3661 

0. 5032 

0, 6353 

0, 7542 ' 

0.8553 

0. 9372 

1 

1. 0453 

1. 0750 

-1. 1 

25 

0. 2333 

0. 3590 

0,4969 

0,6370 i 

0.7709 

0, 8930 

1 

1.0904 

1. 1639 

-1. 8 

27 

0. 1433 

0.2491 

0, 3809 

0, 5306 

0. 6889 

0. 8476 

1 : 

1. 1413 

1, 2684 

-2. 6 

29 

0. 0849 : 

0. 1681 

0, 2862 

0.4359 

0. 6104 

0. 8013 

1 ; 

1. 1988 

1. 3913 

-3. 5 

31 

0. 0484 : 

0. 1104 

0,2108 : 

0, 3532 

0. 5363 

0, 7546 

1 

1. 2635 

1. 5361 

-4. 4 

33 

0, 0266 

0. 0705 

0, 1522 

0. 2823 

0.4672 

0. 7078 

1 

1. 3363 

1. 7069 

-5. 3 

35 

0. 0141 

0. 0438 

0. 1077 

0. 2226 

0.4036 

0.6613 

1 

1.4182 

1. 9091 

-6. 3 

37 

0. 0072 

0. 0265 

0. 0747 

0. 1730 

0. 3457 

0. 6154 

1 

1. 5103 

2. 1491 

-7. 4 

39 

0. 0036 

0, 0156 

0, 0507 

0, 1327 

0. 2936 

0. 5704 

1 

1. 6139 

2. 4351 

-8. 6 

41 

0. 0017 

0. 0089 

0. 0338 

0. 1004 

0. 2472 

0, 5267 

1 

1 . 7306 

2. 7771 

-9. 8 

43 

0. 0008 

0. 0050 

0. 0221 

0. 0749 

0.2064 

0. 4844 

1 

1. 8622 

3. 1877 

-11 



180 200 220 240 260 280 300 320 340 360 


Temperature T, °K 

Figure 3. Temperature dependence the line strengths S? (300). Numerical values are listed 
in Table 4. ^ 




Table 5. Frequency Shift Ti and Intensity Functions f(M’ , N) for Zee man -Components 

of Microwave Spectrum Lines (Eqs. 14, 15), 


Initial State; N, J, |M| < J Final State: U, J', |M'| < J' Selection Rule: AM = M' - M 


Zeeman 

Transitions 

Frequency 
N - Line 

Shift n (m,n) 
N - Line 

Intensity Functi 
(Matrix Eleme 

- Line 

ons f(M' ,n) 
nts) 

K - Line 


M (1 - N) + 1 

M(N + 2 ) + 1 

3 N(N + M' + 1 ) (N + M' + 2 ) 

3 (N + 1 ) (N + M') (N + M' - 1 ) 

(AM = 1 ) 

W(N + 1 ) 

N(W + 1) 

4 (N + 1)^ (2N + 1) 

4 K^ ( 2 N + 1 ) 

(AM = 0 ) 

Mfl - W) 

MfN + 2) 

3 N[(N + if’ ~ M'^] 

3 (N + 1 ) (N^ - M'^) 

n(n + l) 

N(N + Ij 

(H + if (2N + 1) 

( 2 N + 1 ) 

(AM = - 1 ) 

M (1 - N) - 1 

M(N + 2 ) - 1 

3 N(H - M' + 1 ) (W - M' + 2 ) 

3 (N + 1 ) (N - M') (N - M’ - 1 ) 

N(N + 1 ) 

N(N + 1) 

4 (N + 1)^ (2N + 1) 

4 n^ ( 2 N + 1 ) 


After Hill and Gordy ( 195 ^) 

After Lenoir (I968) 


Example : 



n(M,N) 


f(M 

-,N) 


f(N) 


M 


n 

(T 

M* 

0- + 

it 


(Eq, 8 ) 

Line 

+ i 

0.5 

0 

0.5 

+ I 

0,75 

0.75 

0,125 



0 

0.5 

0 

0.5 

0 

0.375 

1 

0.375 

2.5 

J = ^ 1 

- 1 

0.5 

0 

0.5 

- 1 

0.125 

0.75 

0.75 


1 Line 





+ 1 

1 

0 

0 



0 

0.5 

0 

0.5 

0 

0 

2 

0 

2 

J = 0 — J' = 1 





- 1 

0 

0 

1 


































2.4 Line Shapes of the Atmospheric Microwave Spectrum 

The first complete theory of the atmospheric 0^~MS was given 
by Van Vleck (1947). As specific applications based on absorption 
(emission)propertie s of the O^-MS evolved, the need for more detailed 
spectroscopic information was developed. Especially the works of 
Meeks and Lilley (1963), Lenoir (1968), Wilheit (1969), Waters (1970), 
and Westwater (1970) expanded and refined the theory of the atmospheric 
O^-MS for purposes of remote sensing atmospheric temperature 
structure s . 

The molecular theory of the pre ssure -broadened O^-MS was 
advanced by Dillon and Godfrey ( 1969, 1972) and Mingelgrin (1972). 

Both treatments are based on an inter molecular potential model for 
0 -0 , O -N , O -Ar interactions, and scattering calculations of the 

MW W W W 

trajectories for rotational and translational motions. Dillon-Godfrey 
treat the rotational trajectories quantum -mechanically up to N = 7. 
Mingelgrin solves both trajectories classically up to N=23, however, 
the approximate nature of this assumption is expected to yield width 
parameters which are slightly larger than the true ones (R. Gordon, 
private communication, 1973). 

Measured high-pressure (3 to 40 ktorr) absorption rates of pure 
and O^-Ar, O^-N^ mixtures over the 48 to 81 GHz range (Mingelgrin 
et aL , 1972) are reproduced by the theoretical calculations. Implica- 
tions of both theories to the atmospheric O -MS are discussed in 
Section 2.4. 1 . 

In treating the intensity distribution (line shapes F' , F'^) of the 
O^-MS it is assumed that the microwave power flux density is sufficien- 
tly low that nonlinear saturation effects are avoided. As a result, the 
integrated area of each absorption line is proportional to the number of 

^1U 

possible transitions (numbers of molecules) or, on a macroscopic 
‘The number of molecules per cm^ equals 9.662- 10^® p/T (ideal gas law). 


j 
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scale, to pressure. The line strengths (Sect. 2. 3. 2) have been 
defined as scale factors to the area TT • n|^ • y (Eqs. 18, 19), -which is 
that of a Lorentzian shape (Fig. 4). While the theory of spectral 
properties relating to the molecule is -well understood, this is not 
the case for the intensity distribution, especially when considering the 
full range of conditions specified by equation (3). The atmospheric 
O^-MS for pressures above 1 torr is broadened by binary collisions 
(Sect. 2.4. 1); for pressure below 1 torr, the gradual transition from 
collision-to Doppler-broadening is described by a Voigt profile 
(Sect. 2. 4. 3), 

2.4. 1 Pressure Broadening of the Microwave Spectrum 

In the simple case of a single, pressure -broadened line, theory 
predicts the Lorentzian line shape to be valid when y -<< (e. g. , Dillon, 

1969). The frequency, normalized to a multiple of the width y , is 


Z = (Vo - V )/y 

[1] . 

(16) 

yielding the shape factors in the following form: 


Extinction 

Dispersion 


F" = 1/y (1 +z®) 

F' = z/y(l +z®). [1/Hz] 


The approximations for the wing regions (z > > 1) are 

(17) 

F^^^y/(v - V)^ and 

o 

F» ^ 1/ (V^ - V) . 


The peak extinction, ' 

The peak dispersion 


for z = 0 (v = V ) is 
o 

for z = l{v = V ±y)is 
o 


n«» = Sly 
o 

An = S/ 2y 
o 

(17a) 

This yields the two line profiles 

in the form 



n" = n” / (1+z®) and An = An 2z/(l+z'®). [l] (18) 

o o 
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Figure 4 gives a graphical presentation of the equations above. The 
linewidth y can be expressed for atmospheric air by (Liebe, 1969a) 

r=y°(300) im (300/T)'^.p + m (300/T)'^- p }. [MHz] (19) 

Table 6 shows the large body of width parameters y''(300), which exists 

for pure oxygen - 1)- The other quantities of equation (19) are 

and m , the broadening efficiencies for dry air and water vapor, and 
w 

u and v, the associated temperature dependences. The arrangement 
of y°(300) versus the quantum number N in Figure 5 seems to support 
the theoretically predicted N-dependence (Dillon and Godfrey, 197Z; 
Mingelgrin, 1972), although the experimental data (some are claimed to 
be accurate to ± 5 percent) are spread over a large range (± 20 percent). 
The uncertainties of y are one of the largest sources of error in a 
reliable description of the atmospheric O^-MS (Liebe, 1969b). 

The isolated line is a good candidate for quantitative experi- 
mental work on the O -MS. From measured values of An and/ or n^^ 
one needs to determine the quantities (Eqs. 18 and 19) 

U - S°(300). ,3r(T) - , m, , m , u, v . 

Profiles of n'^ and An close to a line center { | - V | ^10 MHz) and at 

low pressures (<1.5 torr) for different values of H can establish the 
validity of equations (14) and (15). O^-MS measurements are in pro- 
gress in our group using dispersion and absorption pre s sure -scanning 
spectroscopy (Liebe et al. , 1973). 

4 

So far, we have measured the 9 line and obtained the following 
preliminary results (Liebe and Welch, 1972): (20) 

4 

(a) V^{9 )-the value is given in Table 2. A pre s sure -induced 

shift of was not detectable (<± 10 KHz/ torr for 50 torr) 
which agrees with predictions (Dillon and Godfrey, 1972). 
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(b) S"'(300)./9'(T) - the agreement with calculated values (Tables 3 
and 4) was within ±5 percent between 250 and 325 More 
careful work at 300 yielded agreement within ± 1 percent. 

(c) F* , F^^— the shape factors were verified as being 

Lorentzian, The pertinent width parameters were 

y^(300) = 1.81 [MHz/ torr ] ± 2 percent, m_=0.93, u = 0.87, 

d 

m =1. 25, and v^ 1 (assumed), 
w 

(d) The experimental An^ versus | - V | profile at H = 0. 53 

[gauss] was in good agreement with '^Zeeman'* calculations. 
The responses of An and n" at H = 0. 53 [gauss] were undis- 
tinguishable from those measured for forp>1.5 

torr and - v|> 10 MHz. 

(e) The comparison between air and oxygen dispersion yielded 

[An^ (Air)/ An^ ^ ^1^ ^d Eqs. 18, 19)- 

For atmospheric O^-MS calculations, we have arranged the rele- 
vant spectroscopic parameters in Table 7 by order of increasing fre- 
quency. The theoretical width values by Dillon-Godfrey (1972) are 
listed, but they have been scaled up by a factor 1. 14 to match our 
experimental value. A computer summation of contributions 

(Eq. 4) from 44 lines of the O^-MS for six different pressures revealed 
the spectra shown in figures 6 through 8. The calculation used the 

S°(300) values of Table 7 and y “ =1.6 MHz/ torr. Pre ssures chosen 

1 

are representative for h^ 30 to 0 km. We conclude from these results 
that the atmospheric O^-MS needs to be discussed separately for three 
pressure ranges: 

I - Isolated Liine Spectrum for p< 10 torr (Fig, 6). 

II “ Continuum Spectrum for p > 300 torr (Fig. 8), 
where the envelope is determined by all lines. 

I/II - Mixture of I and II (Fig. 7). 
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In range II, the addition of individual line responses (Beer- 
Lampert law) may not be a valid procedure since collisions do alter the 
energy levels to an extent that they can interfere with each other. 
Mingelgrin (1972) included overlap effects in his calculations, which are 
compared with the pre s sure -linear width model (Eq. 19) in Figure 9. 

A substantial narrowing can be noticed above 100 torr. For standard 
conditions (760 torr, 298 "K) the difference in peak intensity at 60 GHz 
is as high as +28 percent. Overlap effects between N and N lines 
^have been computed by Dillon (1969) and found to be insignificant for 
standard conditions. These contradicting predictions need to be 
clarified through analysis of reliable experimental continuum data. 

A necessary condition for overlap in the coalescence of extinc- 
tion line shapes (Eq. 17) under the influence of pre s sure -broadening. 
Table 7 reveals 8 lines where the separation 0 of the centers is only 
about 1 / 5 of the normal spacings. These doublets overlap with 
increasing pressure first. A second condition comes from pressure - 
broadening theory, which predicts different sensitivities of the energy 
levels to mutual interferences. Approximately, one can treat the 
individual line widths as being no longer proportional to pressure and 
formulate 

y^(range II) ^ y.” p - 6 ( 0 , p, N) . ( 21 ) 

For example, the 760 torr attenuation profile by Mingelgrin 
(Fig. 9) can be fitted approximately by assuming for the doublets a 
line width of one -half of the isolated width. Thus the doublets may be 
responsible for most of the overlap effects of the atmospheric continuum 
O^-MS. Considering only overlap factors 6 for these line pairs would 
simplify computations of the continuum O spectrum. 

An empirical width model was introduced by Meeks and Lilley 
(1963) to reconcile atmospheric attenuation rates between the line 
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spectrum (I) and the continuum spectrum (II), Zenith attenuation data 
(see Table 10) have been fitted with this model whereby over the range 
of the continuum spectrum, >' was reduced to 1/ Z (Carter et al. , 1968) 
or even to 1/ 3 (Reber, 197Z). This will be treated further in Sections 
3. 1 and 2. 

2.4.2 Nonre sonant -Spectrum 

In a classical picture the molecular dipole is unable to follow 
completely the oscillating field. Thus, after every collision there is a 
preferred reorientation relative to the field. Formally, collisions in- 
duce transitions between identical energy levels under the selection rule, 
AJ = 0 and AN = 0. As a consequence, increasing pressure (p> 100 torr) 
adds to the O^-MS a non-re sonant contribution which is centered at 
V = 0 Hz and can be described by a Debye -type polarization (VanVleck, 
1947; Zhevakin and Naumov, 1967) which is 

n" = S F" and An = S F' , (22) 

n o o n o o ' ' 

The Debye shape factors are 

F" = 2 (v®+y®) and F' = jy ® / (v® + y^) J - 1 . (23) 

The dispersion spectrum at 760 torr was calculated (v= 0 to 80 GHz; 

N ^ 25)by Mingelgrin (1972), and from the result, which is shown in 
Figure 22, we can deduce approximate values of strength and width 

2. 10'"° (300/T)^.ri.p , [l/torr] , (24) 

y^ 1. 8 (300/ T) . p , [MHz/ torr] . 

The example for standard condtions at sea level (760 torr, 300 °K) 
shows that at microwave frequencies the nonre sonant part is small when 


28 



compared with the much larger resonance part and it can very well be 
neglected in the remainder of this report: 


V [GHz] 

1.37(y^) 

40 

90 

140 

a [dB/km] 
n 

. 021 

. 041 

. 041 

. 042 

A 0 [ rad/ km] 

n 

-. 001 

-. 070 

157 

-. 244 


2.4. 3 Transition to Doppler Broadening (Voigt Profile) 

As the pressure decreases with altitude (Eq. 3), we have to 

consider another broadening mechanism - the Doppler effect. The 

transition frequency of a non-colliding molecule moving with velocity 

c in respect to the radiation field is v [ 1 + {c / c)l . The velocity 
m o m ^ 

distribution of gas molecules in thermal equilibrium is Maxwellian. 

The resulting Gaussian intensity distribution has a linewidth which 
depends on temperature. The pre ssure -independent width makes the 
peak intensity proportional to pressure, thus forcing the atmospheric 
line spectrum to vanish. The equation for (see Table 8) yields 
the following values for 


T["K] 

y^[kHz] 


At V - 68. 49 
o 

and at 49. 96 [GHz] 

180 

59. 1 

42. 5 

290 

74. 9 

53. 9 


using 

y^^ = 6. 34 • 10“® • •V’t'’ . [Hz] (25) 

The Gaussian shape of a Doppler -broadened line is for extinc- 
tion (F") a little broader at the center but decreases much more rapidly 
in the wings when compared to a Lorentzian. A normalization factor, 


29 




Vn In 2 is introduced for the Gaussian (area = n" y- '/n/ln2 ) in order 

o 

to match the Lorentzian (area = n^^.y.TT) integrated strength (Table 8). 

A convoluation of the two shape functions leads to a Voigt profile 
whereby it is assumed that the free motion of molecules is unaffected 
by collisions (’*pre ssure and ’Temperature'’ are statistically independ- 
ent) which is not exactly true. The reduced mathematical expressions 
for the Voigt function are given in Table 8. The Voigt parameters 

y=y°p/y^ and x = Av/yj^ (26) 

govern the extinction profile u(x, y) = n”/ n|^ and the dispersion profile 
v(x, y) = An/ An^ as shown in Figure 10. For y< 0. 1 the shape is 
practically Gaussian, and for y> 10 the shape is Lorentzian, As the 
pressure goes to zero (y<0. 1), the rate for the Lorentzian values nj^ and 
An to approach zero becomes 0. 678 S°p/y^ for the max. extinction 
n^ and 0. 610 S®p/ for the max. dispersion (see Figs. 10 and 11). 

Figure 11 also gives the corresponding frequency loci of An^ (p) and 
y(p). In principle, the Voigt profile has to be applied to each Zeeman 
component in the presence of a magnetic field H. 


2. 5 Attenuation and Dispersion Due to Water Vapor 

Two rotational lines of ^ 2 ^ centered on either side of the 
40 to 140 GHz band. The low-frequency wing of the complete H^O 
rotational spectrum and additional polarization mechanisms also 
contribute to atmospheric millimeter wave transmissivity. We can 
separate the resonance contributions from the remainder (q) and define: 


2 2 

n”(v)=:E S. F;' An = i: S F: +q’ 

w i^i^i w w 1 1 w 


(27) 
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The resonant contributions by the two lines adjacent to the frc ^ 

quency band under consideration can be evaluated using the following 
spectroscopic parameters (L.iebe, 1969a) 


i 

1 1 

1 

V 

o 

[GHz] 


22. 236 15 



183, 310 



S 

[Hz] 

13, 9( 

3 6 

r300 \ r . . , 

r, 300\1 

,,,/300' 


300\1 


■ expj^2.14| 

[■—)]'>» 


j ■ exp 0 . 6 5 3f 1 - 

T /J 

•Pw 

y 

[MHz/ torr ] 

18. o' 

[o.209(^) • p+( 

300\ I 

“)PwJ 

19[o,2l( 

300\'^'^ /300\ 

— ) ■ pn~) 

■Pw] 



The Gross shape was found to give a better fit of experimental data in 
the wings of strongly foreign-gas*-broadened lines, 

F> = 2v (V® - v^)/ [(v® - + (2vy)^] , [1/Hz] 

Lx O O O 

(28) 

F"!. =4v®y/[(v^-v^)®+(2vy)®] . [1/Hz] 

Lx o 

An example for p - 10 torr under standard conditions, p = 760 torr and 
w ^ 

and 300 ®K, yields the following resonance contributions (Eqs. 2, 27, 28) 


v[GHz] 

(V ) 

22^ 40 

60 

80 

100 

120 

140 

183^^0^ 

a [dB / km ] 
w 

. 194 

.016 

. 009 

. 010 

.017 

. 038 

. 106 

33. 8 

A0 [rad / km] 
w 

. 001 

. 026 

. 047 

. 071 

. 103 

, 153 

, 245 

. 017 


These contributions are by far not enough to account for milli- 
meter wave attenuation. Laboratory measurements at two fre- 

quencies produced results shown in Figure 12. Experimental uncer- 
tainties are on the order of 10 percent. A tentative fit of the data was 
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accomplished by (Liebe et al, , 1973) 


(^ 9 ) 


2 4 2.6 10 


where the frequency v is in gigahertz. 

The term linear with p indicates roughly the contribution 
expected from the residual water vapor attenuation due to the low fre- 
quency wing of the rotational spectrum. VanVleck (1947) estimated for 
this term the following amount (Liebe, 1969a) 

a 1. 9 • lO"^* V ^ • p / T^' ^ , [dB/km] , 

w 

_3 

which yields at 6l GHz and 300*K the value 4. 5 • 10 [dB/km torr]. 

2 

The p^ term, however, points to pre s sure -induced polarization 
effects such as dimerization. Noteworthy are the temperature depen- 
dence. An Arrhenius plot, loga versus T \ yields a "binding energy" 

of 4. 5 [kcal/mole] ± 10%. The p^ term gives rise to abnormal high 

2 

attenuation when saturation (pg) is approached. A similar p response 

was measured for the dispersion q' = n' (v.) - n' (v,/ 2) for v, at 

61.2 GHz (T = 280 and 325 *K). The sketchy experimental evidence of 

unaccounted, pressure -induced polarization in water vapor certainly 

warrants a more carefial study of a (q" ) and q' over a wider rang of 

w w 

millimeter wave frequencies and close to saturation pressures p (T). 
Such experiments can shed light on the fact that the theory of the rota- 
tional spectrum does yield atmospheric attenuation values in the 

wings of the two microwave spectral lines which are too low when 
compared with measured data. 
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Normalized LORENTZIAN Line Shape 




Table 6, 


Reported width parameters y^(300) and their temperature dependency u for self- 
pressure broadened lines of the Microwave Spectrum. 


4 ^ 


N± 


(c) 

ijd) 


yO{30i 

E 

(f) 

3) in 
xperir 
(b) 

MHz/tc 

nent 

(i) 

>rr 

(k) 

U) 

(b) 

, (c) 

u 

(d) 

(e) 

(k) 

(f) 

1 + 



CO 

00 

a- 

r- 



1.90^®^ 


2. 20 

1.97 

1, 96 

.80 


o 

i 

ID 

00 

d 

.90 

“ 

3f 



1. 86* 

1. 96 
1.96 

2.23 

2. 07 

1.71 

1.92 

.81 


. 85 

2 

5+ 



1. 72^ 
1. 98* 

1. 56 
1.60 

1.96* 

1.99* 

1. 80 

1.86 

. 79 


. 88 

2 

2 

7 + 



1. 73* 
1.53* 

1. 68 
1. 70 

1.92* 

1.82* 

2. 01 

2. 05 

.79 


. 84 

2 

2 

9+ 

1.95 ( 1.59) 

1.81 

1. 78* 
1.48* 

1.42 

1.64 

1.93* 
2. 00* 

1. 94 

1. 97 


. 85 

.71 

2 

2 

11 + 

1.95 (1.59) 


1.56* 

1.60 

1.97 


1.97 




2 

13 + 

1.91 (1.59) 


Reduced from 

310 to 300’’K 

1 1 1 

1.54 

1.86* 


1.87 





15 + 

HHIHH 


— 

1.99 


1. 77 





17 + 



1. 50 



1. 76 





19 + 

1.58 (1.59) 

i 


1.92 


1 . 

1.62 





21 + 

— If 

1 — ■ 


1 r* ^ 



1. 26 





23 + 

BBI 


X 

X D percent tiiio 

1 L_ 

l."2'^” 

1.29 






References: (a) Mingelgrin (1972) (c) Liebe Welch (1972) (g) Schulze Tolbert (1963) 
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Rotational Quantum Number N 

Figure 5. Reported line width parameter yt^OO) for the self -pressure broadened 
Microwave Spectrum (see Table 6) versus the quantum number N 
for end-over -end rotational angular momentum. 
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Spectroscopic parameters of the O 
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microwave spectrum 


used for computations of this report. 
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Figure 6-8. Computer plots of the pres sure -broadened Microwave Spectrum for six 
different pressures (10 to 760 torr) corresponding to altitudes between 
and 0 km. The model of linear addition of individual Lforentzians was used 
{T=300^K, v! = 1. 6 MHz/torr). 

For analysis of atmospheric conditions see Figures 14 through 17. 
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Figure 9. Attenuation of oxygen between 54 and 66 GHz at four dif- 
ferent pressures. 

Solid lines: Linear addition of Lorentzians (see Figs. 6-8) 
using Mingelgrin* s width parameters (Table 6). 
Dashed lines: Scattering calculation by Mingelgrin (1972), 
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Table 8. Shap 


Line Profile Broadening Mechanism 


LORENTZIAN Pressure p 

(Collisions) 


oo 


(See Fig, 4) 

VOIGT Convolution between 

Gaussian and 
Lorentzian 


(See Fig. 10, 
y 0) 

GAUSSIAN Doppler Effect due to T 

(Velocity distribution). 

Doppler width, 


Yd= (v^/ c) V a kT/ M 


M - molecular mass 


(a) 


Miller and Gordon (1931) 





1 + z / ^ 


^Frequency’* z = (v^-v)/Y 


“(==.y)=yj' 


"Frequency” x= 
"Pressure" y ~ 

I 


4VI^ r , , =1 = 

F =— exp[-ln2x J J exp (x’) dx' 

» D o 


VnlnZ / 1 

p ~ — exp (-InZx ) 

^D 


Dawson’ s Function 





4757 



Figure 10^ Normalized Voigt profiles for dispersion v(x,y) and 

extinction u(x,y) (Faddeyea and Terent'ev, 1961). The 
mathematical functions are given in Table 8. 
"Frequency" x= (vo- * solid-line curves. 

"Pressure" Y = * dashed-line curves. 
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Frequency Ai/ , Units of 



Pressure, Units of p 


Figure 11. Frequency loci of the maximum dispersion | An | 
and line half -width (y) for a Voigt profile. Also 
shown is the normalized amount of | An^ | . The 
ordinate is in normalized pressure units (y). 
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Attenuation 0! ,dB/km 



Figure 12. Pure water vapor attenuation as a function of the vapor pressure for two different 
temperatures at 61. 16 GHz, and at 30. 58 GHz. 




3. COMPUTER ANALYSIS OF ATMOSPHERIC PATH MODELS 
3. 1 Homogeneous Path Transmissivity at Various Altitudes 

Before applying theoretical results of Section 2 to actual atmos- 
pheric conditions, it seems appropriate to establish their validity. A 
serious test of the theory is possible by comparison with experimental 
results of sufficient accuracy which in turn can only be obtained under 
controlled laboratory conditions. 

The Line Spectrum 

One experiment yielded quantitative attenuation data for eight 
isolated lines around 60 GHz (Stafford and Tolbert, 1963) and the single 
line at 119 GHz (Schulze and Tolbert, 1963)* The conversion of peak 
attenuation 0,^ and linewidth y° to a strength parameter (Eq. 11) is 
accomplished by means of (Eqs. 2, 18) 

S° = 5493 a y"/v [Hz/torr] (30) 

o o 


where is in megahertz (from Table 7). This is illustrated by the 
tabulation below 
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The difference between S® and is for five lines larger than 

the combined experimental uncertainties. The fit improves when the 
linewidth parameters of Table 7 are applied. The discrepancies are 
probably due to insufficient control of experimental conditions, incor- 
rect assumptions in reducing the data to 310°K (u = w = 2 !), and baseline 
(a = 0) problems. This, together with our experimental results for the 
9^ line (Eq. 20), leads us to believe that the theoretical line strengths 
are valid as given by equations (7), (11), (12), and Table 7. 

The other important parameter of the isolated line spectrum 
is the correct width. Ab initio calculations of y“(300) by Mingelgrin 
(1972) and Dillon -Godfrey (1972) are sensitive to the intermolecular 
potential parameters which are obtained from data unrelated to the 
O^-MS. Results from both theories differ by about 20 percent as do 
the experimental values (Table 6, Fig. 5). Confidence in the measured 
y®(9 ) parameter (Eq. 20) and support by other experiments of an 
N-dependence for the width (Fig. 5), made us use the 1. 81 MHz/ torr 
value to "calibrate" the Dillon -Godfrey widths. 

The atmospheric transfer function versus frequency above an 
altitude of about 25 kilometers resembles a string of individual 0,-MS 
lines. The analysis of the line spectrum seems to be in fairly good 
shape. Uncertainties are connected with the width parameters 

19) for other lines than the 9"*^- Further 
experimental work with the dispersion spectrometer (Liebe et al. , 1973) 
will hopefully provide the missing information. A guide to treatments 
of the line spectrum in this report with respect to practical applications 
can be found in the overview on the next page 
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Guide to Treatments of the Atmospheric Oxygen Line Spectrum (h > 25 km) 


Altitude 

Line Broadening 
Line Shape 

h 

> 25 km 
Pressure p 
Lorentzian 

> 40 km 

Magnetic Field H 
Lorentzian - Voigt 

> 60 km 
Temperature 
Voigt 

T 



Eq. 

Table 

Fig. 

Eq. 
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Fig. 
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1.2 
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10 
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5, 13 

19 

7 
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The complete O^-MS aad detailed insight into its transfer pro- 
perties require a good deal of theory. Therefore, it is important to 
test theoretical results by as many reliable experiments as possible. 

Interesting properties of isolated lines are their markedly 
different temperature dependences. The peak attenuation varies 
approximately as (Eqs. 2, 13, and 18) 

a - (T /T)^“^ . (31) 

o o 

Looking at Table 4 and assuming u = 0. 9, we find that of the two 
N= 15 lines is almost independent of temperature and thus a good 
candidate as an atmospheric pressure sensor (from a satellite). One 
can also pick lines with inversely linear (N= 11) and square -law (N= 3) 
dependences as well as linear (N= 15), square -law (N=23), and higher 
order temperature dependences. The high rotational lines are very 
sensitive functions of temperature (Waters, 1973). For example, a 
1 percent ( + 1. 8 °K) raise at 280 produce s for N = 27 line s a 3 . 5 per - 

cent increase in , and at 220 the increase is even 5. 1 percent. 

The Continuum Spectrum 

One early experiment (Artman and Gordon, 1954) provided 
attenuation data (i 10 percent) at 25 frequencies between 52 and 6l GHz 
under controlled laboratory conditions (p = 760, 380, 190 torr; 

T = 300‘^K). Several workers (e.g., Meeks and Lilley, 1963; 
Westwater, 1970) used these data, different line shape functions, and 
theoretical line strengths to deduce the pre s sure -broadening line 
parameters. It was found that the isolated width parameter (assumed 
equal for all lines) needs to be reduced by 30 to 50 percent (see Eq. 12) 
to fit more or less the continuum spectrum, and in all cases it was 
stated clearly that future laboratory measurements are needed to 
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remove the encountered uncertainties. Since such experiments have not 
been reported up to now, we have had to rely upon measurements made 
in the atmosphere. 

Results of field measurements contain information not only 
about the O^-MS but also about the structure of the atmosphere. Many 
of the uncertainties in present knowledge of the transfer fvmction cannot 
be distinguished from uncertainties in pressure and temperature distri- 
butions which makes field data an unsatisfactory alternative to laboratory 
measurements. 

Once the O^-MS is completely understood it is an easy task to use 
spectroscopic properties as a probe of (dry) atmospheric conditions. Right 
now, however, we have to content with an empirical line width model 
introduced by Meeks and Lilley (1963) which determines the individual 
width from (y. applied to all lines except the 1 ) 

y^ (p) ^ G(h) (p/ 760) (300/ T)'^ and (32) 

J y I (p) + 

Carter et al. (1968) and Reber (1972) measured oxygen attenuation of the 
atmosphere by the zenith angle technique using the sun as a source. 

Their ground level results are given in Table 10. Carter et al. accom- 
plished a ± 8. 8 percent fit (squared deviations) to their data which 
include airplane measurements (taken between h = 2. 6 and 12. 2 km) by 
the following scheme for G(h); 
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G(h) 
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[torr ] 

y.“(300) 

[MHz/torr] 

0 to 8 

640 

0. 85 

>180 

0. 84 

8 to 25 

640 to 1357, linear with h 

M 

- 

- 

> 25 

1357 

rr 

<20 

1. 79 
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With these values, equation (3Z) becomes the ’*CMR-Model'^ which is 
used for the analysis described in the remainder of the report. 

Another approximation, equation (33), is introduced to take care 
of the transition to Doppler -broadening (Section Z. 4. 3). The Lorentzian 
shape is assumed for the vanishing O^^MS, and the transition to a 
Gaussian shape (Voigt profile) is not done. At h = 50 km this approxi- 
mation was found to be indistinguishable over the strongest 90 percent 
of a line from the more correct Voigt profile (Waters, 1970). Equation 
(33) is sufficient in the context of analyzing cumulative effects through 
the atmosphere (Sect. 3, Z), but for a calculation of the line profile over 
a homogeneous path above h = 60 km the Voigt profile has to be applied. 

3. 1. 1 Analysis of Attenuation and Dispersion for Dry Air Between 
49 and 7Z GHz, 

The analysis includes all lines with rotational quantum numbers 

N= 1 to 43 (Table 7), which omits lines with strengths less than about 
-5 

Z • 10 of that of the strongest line. This cutoff is sufficient even for 
computations of transmission over long paths. The assumptions made 
for the computer analysis of homogeneous attenuation and dispersion 
rates (Eq. Z) were: 

(a) The U. S. Standard Atmosphere (196Z) models the distribu- 
tion of pressure p and temperature T versus altitude h. 

(b) The atmosphere is dry, p = 0. 

w 

(c) The sum of 44 lines (Eq. 4) is the only source of attenua- 
tion and dispersion over the frequency band considered. 

(d) The intensity distribution of each line is a Lorentzian (Eq. 
(Eq. 17), and the theoretical strength values (Table 7, 

Eqs. 7, 11, IZ) are valid. 
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(e) The individual widths are the same for each line. They 
are given by the *^CMR-Model*^ (Eqs* 3Z, 33) whose 
approximate (empirical) nature provides at best analytical 
estimates of transfer properties. The altitude dependence 
of the width is shown in Figure 13. 

Results of the calculation are given in Figures 14 through 19. Only the 
graphical presentation provides an overview of the structure of the 
transfer function at different altitudes. Numerical results (computer 
print -outs) for any given set of conditions are available upon request. 

At h=0 km, figure 14 shows a continuum spectrum (II). The 
envelope is not a sensitive function of the width parameter: 

(1) y° = 0,876 MHz/ torr plus arbitrary variations of ±25 percent (CMR- 
Model); (2) = 1. 27 MHz/ torr, a value obtained by Westwater from 

fitting Artman and Gordon* s (1954) data. Atmospheric attenuation for 
sea level conditions was reported at 21 frequencies between 40 and 120 
GHz (Thompson, et al, , 1972, Fig. 3). 

At h = 10 and 20 km, figures 15 and 16 show a mixture of 
continuum (II) and line (I) spectra. The influence of the width becomes 
more pronounced as the lines gain identity. At h = 30 km (Fig, 17), 
we have a pre ssure -broadened line spectrum. 

At h = 50 km, figures 18 and 19 show the effects of Zeeman 
splitting of two individual (9^ and 27 ) lines. Figure 18 gives changes 
in. transmissivity for the case of linearly polarized radiation. Each of 
the 19 TT - or 38 Q -components (Eq, 14) is treated as a pre ssure - 
broadened line of Lorentzian shape (Eq. 17) with a width determined 
by equations (32) and (33), The line strengths are calculated using 
equation (15). The contributions from each Zeeman component are 
summed up over frequency (Eq, 4), The 9"^ line widens under the 
influence of H = 0, 53 gauss while the peak attenuation is reduced by 
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ZO to 36 percent, respectively, for the TT - and a -envelope s. The 
a -envelope ’splits for circular -polarized radiation ^spending on the 
orientation., . The^^two halves are shown in the, exanipl^_ the Z7 line 
in figure- 19 .:- bXhe ind>yiduarZeernan cpr^ ^re,n(^t resolved in 

both cases. A consequence of the ha 40 km 

the atmosphere-.bec.omes a medium with polarization, dependent aniso- 
tropic transfer^ properties. , , . ^ ■ 

. „The evaluation of a practical, case, requires the.,d^te.rrninatiqn.of 
the type-of polarjizatiqn for the propagating. radiation. | The ne:^rt. step is 
to define the orientation between the direction of, H and,,the.. miqrowave 
field component_^J^,,. Lenoir/ s (196§) transfer function in tensor form 
does this in general form for geomagnetic coordinates, Cr^oqm (1971) 
disGusse|d 24, posisible "pure,'! constellations between H and for, the , 
case, of- emission from the 1 line. whe re; ^ takes the .rqle of the radio- 
meter antenna polarization. 

The simple Zeeman pattern of fhe 1 and 1"*^ lin.e.s ease 5 . 
considerably the i-nterpretatipn of emission (atten.ua;tion) measurements 
in. the vicinity, of the.se. lines. . Such emissions can be picked up, from 
satellite jborur;riadi ome ter s an, d .then .cohy^rted t^o me sospheriCj.ternpera- 
ture structures and magnetip fie.l^ (^) data , 00 ;, a. global- h,asis (Lenoir, . 
1968.;-,Crpom,i^. 1971,)..- At the se heights the L"*^jline is-woll-isole^ted from 
itS;:l!^n^:'5at®^%t^•I T)^® jP?ak atte>nUMion,of bo,thjline,s- (Eq.s., ,2,„ 7,, 1.8,;. , 
andiTable 7): have for-;T 300 v-K the, yalu^^^ v.ii. .'rlc^£,^ o. 


"'a - ( 1 '+)-'=' : 35 Land a • (i“)> = 1; 3 &-^'tdB/-km 0 -..o - v; l i 
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‘ Crd'bm-’( l'9flV'g^^'e'^ li i-0‘for^ M8‘*K‘* ^ a valne"1x^’( i-^) ’=-^3P8' dWVih' 
whiG)j. is WfTong-sixxc^G-tiie;' o\^;rlap-,red^ce^ (-^q.v j3Z) was -jii-sed;. At 

288 we calculate 1. 54 dB/km (Eqs. 2, 11, 12, 18, 19, and Table 7). 
This ^err6y i‘e(i'uce's"the zenith^ '^tteriuatib whibh shoiild clfeWex^td' 
Aq(T T-).^ i lp- dB (p. 7* 1-^5* 100,km);djistG^^ of. 26^ dB jdiS^ givp^n 
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Figure 13. Linewidth y. (same for alHines of -MS) versus altitude h and pressure p as 
obtained from the CMR -model (Carter et al. , 1968). 
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Figure 14* Horizontal (homogeneous) transmissivity at sea level, h =0 km 
Variations are shown due to different linewidth values : 

(I) 7 . =666 MHz ± 25%, 0 Y. =968 MHz. 










Figure 16. Horizontal (homogeneous) transmissivity at h-ZO km 
a.U) inA 01 =(%.%. (Th«3]ahe^dth.vflsliheSJlaE-e-l y..':^.8:6'i4'MHz 
-anil ineraiiib oj aul' srtohsr .-' 'I d i ■ ■, ^ 
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Figure 18. Effect of Zeeman -splitting upon the horizontal (homogeneous) 
transmissivity at h = 50 km (U. S. Std. Atm. 62) in the vicinity 
of the 9"^ line of the O -MS for linearly polarized radiation 


with *^1 H (tt) or It H {(T). 
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Figure 19. Effect of Zeeman -splitting upon the horizontal (homogeneous) 
transmissivity at h = 50 km (U. S.Std. Atm. 62) in the vicinity 
of the 27“ line of the O -^MS (after Waters, 1970) for linearly 
polarized radiation with *^_LH(n) or for circularly polarized 
(right-, left-handed) radiation (a"^, Q“). 




The intensities of both lines differ by about 1 to 4, but technology at 
56 GHz (!■*■) is more advanced. The NIMBUS 5 radiometer, for 
example is working well in space at 60 GHz (D. Staelin, private com- 
munication, 1973 ). 

3, 1. Z Phase Dispersion Between 10- and 140- GHz 

While attenuation (or emission) is the subject of all papers 
dealing with applications of the O^-MS, data on phase dispersion are 
not readily available. In general, the complete transfer function 
T (cc , 0) (Eq. 1) is required for an evaluation of atmospheric limitations 
imposed on system performances. The phase term has particular rele- 
vance for broadband communication and for various disciplines of micro- 
wave measurements. For example, phase dispersion limits the useable, 
phase -coherent channel bandwidth of a broadband system (Morgan and 
Ekdahl, 1966); propagation delay and its fluctuations determine the 
accuracy of clock synchronization via satellite (Levine, 1970), the reso- 
lution of long -baseline interferometers (Schaper et al, , 1970), the 
accuracy of radio range measurements (e.g., Thompson, 1968), etc. 

The atmosphere delays radio signals causing an apparent 
increase AL in the geometric path length L. The propagation time 

t = ( L + AL )/c [s] (34) 

is determined for refractive indices n' ^ 1 by the phase velocity c/ n' 
and for n'< 1 by the group velocity cn'. 

The case n'< 1 exists for the O -MS only under special circum- 

Ld 

stances when the negative resonance dispersion An(v) of an isolated line 
(Fig. 4) is larger than the non-dispersive refraction (Section Z, 2). 
This case occurs at frequencies up to several megahertz above a line 
center (Liebe et al. , 1973, Fig. 4), but it will not be considered in the 
following discussion. 


58 



We define 


AL = J (a' - 1) dl = J [N- 10"^+ An(v)] dl [ m] (35) 

L L 

where the integration is taken along the ray path. For example, the 
additional electrical length of a zenith path from sea level to outer 
space is typically AL ^ Z meters (Schaper et al. , 1970), Extreme 
values of zenith path length changes due to O^-MS dispersion are only 
+ 1Z millimeters at 58. 3 GHz and -14 mm at 6Z. 6 GHz. The extreme 
values for AL(v) from different initial altitudes are listed in Table 9- 
The calculation procedure of these results will be described in Section 
3, 2, 

System engineers express dispersion effects of the homogeneous 
atmosphere by a phase change in radians per unit distance; 

A0(v) = (2ttv/c) An(v) [rad/ km] * (36) 

Calculations of A0 over the 49 to 72 GHz range (O^-MS attenuation) have 
been described in Section 3. 1, 1, and various examples for h = 0 to 
50 km display the O -MS dispersion spectrum (Figs. 14 through 18), 

It is evident that the molecular resonance dispersion extends over a 
much larger frequency range. The reason lies in the shape function F* 
whose wing intensity decreases with | v -v| as compared to | 
for F” (Eq. 17 and Figs. 4, 10). 

The calculations are extended to cover the frequency range 10 
to 140 GHz for h ~ 0 and 10 km, and the graphical results are shown in 
figure 20. The dispersion of the 22 GHz line puts a ’^wrinkle” into 

the dry air phase response. At h = 0, we calculate (Eqs. 27, 28, 36): 

_3 

(A0 ) =: ± 1. 2 - 10 rad/ km per torr of water vapor at 19- 2 and 

w max 

25. 3 GHz respectively, while the -MS dispersion over the same fre- 

“3 

quency range is + 25 • 10 rad/ km. 


59 



The dispersion An(v) of the O^-MS was first calculated for h = 0 
by Zhevakin and Naumov (1967). The complete oxygen dispersion 
spectrum including overlap and nonresonant dispersion (Sect. 2.4) was 
recently evaluated by Mingelgrin (1972), and his result at h = 0 km 
for V = 0 to 80 GHz is shown in figure 21. The peak value span (7. 87 
ppm) of translates for air into 2. 1 1 rad/ km, which agrees well with 
the value of (2) in figure 14. Another dispersion calculation for dry air 
at h = 0 km was performed by Sullivan and Richardson (MITRE, 1965). 
They predicted a dispersion of 0. 4 ppm between 45 and 90 GHz while 
our calculated value is 0. 381 ppm. A broadband system operating on 
the ground close to the O^-MS attenuation band, e.g. , in the 45 or 75 
GHz bands, does experience phase dispersion on the order of ±0. 03 
rad/ km per gigahertz. 

Dispersion measurements are difficult in the laboratory because 
the small effect is heavily masked by instrumental instabilities, A few 
measurements have been made in our laboratory using a differential 
microwave refractometer (Thompson and Vetter, 1968) similar to the 
one described by Liebe et al, (1973). The difference of the refractivity 
at Z6 and at (26/ 3) GHz was measured at 610 torr and 300 °K in dry air 
to be 0. 12 ± 0. 03 ppm (M, J. Vetter, private communication, 1968), The 
calculated value for the same conditions is (0, 394 - 0. 261) = 0. 133 ppm. 

In another experiment we measured the difference in refractivity 

between 23, 6 and (23. 6/ 2) GHz (Liebe, unpublished results, 1969). In 

dry air from p = 50 to 600 torr at 300 ®K, the following result was obtained 
-4 

( 1, 25 ± 0. 3)' 10 ' p, while the calculated response is 

(6, 03-4, 55) - 10 ^*p = 1.48* 10 "^*p [ppm/ torr] . 

These two experiments provide a rough (±25 percent) confirma- 
tion of O^-MS dispersion outside the 60 GHz attenuation band where the 
dispersion is nearly independent of the widths y (Eqs. 19, 21). 
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Atmospheric refraction limits the accuracy of such systems as 
©lectronic distance meters and tracking radars. The measurement of 
distance is based on determining the propagation time t (Eq, 34), 
whereas the deviation of the vertical tracking angle is proportional to 
the gradient dn’ / dh. These errors can be reduced if pertinent atmos- 
pheric parameters of the propagation path are known. The dispersion 
of the O^-MS offers a means of evaluating the refractive index n' of dry 
air averaged along an atmospheric path. Systems operating at two 
(or more) phase -coherent frequencies can extract path conditions from 
the difference in phase shifts experienced by the signals when travelling 
the same path (MITRE, 1965; Thompson, 1968; Liebe, 1969b). 

The potential of a differential phase measurement is analyzed 
for signal frequencies in the wings of the O^-MS. Two different sig- 
nal pairs are chosen for short and long operating ranges in the following 
example 
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The phase shifts 9 experienced by each signal over a path length Li 

i , L 

are (Eqs. 35, 36 ) 

= [2nVj (1 + N- 10'^)/c + A 0 (v^)] . L 

and , (37) 

02 = [2 nv^(l+ N-10 )/c - A0{v^)] • L , 

when the signal frequencies are chosen so that |A0 (v )| *=» | A 0 (v ) | 

JL ^ 

(see Fig. ZZ) and N is the refractivity averaged along the path. The 
averaged refractivity can be separated into a 'Mrv^^ and a term, 

N = + N 

d w (38) 

Assuming (a) a straight line -of-sight ray path (no reflections or 
bending), (b) a perfect correlation of time -varying phase changes 
(mainly due to water vapor density fluctuations), and (c) a common fre- 
quency v = (v^ + ^ 2 )^ 2 (^' g* » t>y heterodyning and and then multi- 
plying each IF by the opposite factor X) one can measure at the receiving 
end the difference 

” ®2 “ 0L =|a0j(v) + 1 A 02 (v) 1)> *L [rad] ( 39 ) 

and the sum 

+ 9 ^ 2n (Vj + v^) (1 + N 10~^) L/ c [rad] . (40) 

The O -MS dispersion 0 can be tabulated in terms of dry air density, 

3 

= 464 p/ T [g/ m ]. Approximate pressure and temperature 
sensitivities are listed on the previous page. For the path length L, 
first the uncorrected value can be used with some possible later iter- 
ations. Thus the value found through a measurement of 0 determines 
the dry term refractivity (Eqs. 5, 39) 

= 0. 223 [ppm] . (41) 
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Table 9. Extreme values of zenith phase dispersion AT and electrical 
path length change ALr(v) with associated frequency v and 
attenuation A for different initial altitudes, h^ (CMR-Model). 
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Figure 20. Phase dispersion in dry air due to the microwave spectrum of O for horizontal paths 
at h = 0 (760 torr, 15 ®C) and h = 10 km (199 tor, -50®C) using the linear model (sum 
of Lrorentzians). 
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Figure 22. Phase dispersion A0 and associated attenuation a versus 
temperature T for two different harmonious frequency pairs 
in dry air at sea level using the linear model (sum of 
Lorentzians). 0 - differential dispersion between 5 1 . 0/71. 4 

and 55/66 GHz. * 



The wet term refractivity N originates from the water vapor 

w ^ 

within the path. With equation (5) and p^ = 288. 8p^/ T [g/ m ] one 
obtains 

N ~ 1. 73 • 10^ • ■p / T [ppm]. (42) 

w w 

It is possible to solve for the average water vapor density along the 

path, p using the measured sum of phase shifts (Eqs. 40, 38, 41, 4Z). 
w 

The proposed measurement principle offers the potential to 
follow air and water vapor density fluctuations separately in real time 
as they occur as spatial averages along a line -of -sight path. This 
information is also a wanted correction term (Eqs. 34, 35) for radio 
distance measurements. 

3. 2 Slant Path Transmissivity 

3.2. 1. Calculation Procedure and Critical Comments 

Computations of non -turbulent transfer properties through the 
first 100 kilometers of the atmosphere depend on the exact functional 
form of a(h) and A0(h) [Sect. 2], on the model atmosphere which deter- 
mines p(h) and T(h) [U. S. Stand. Atm. 1962], on the mixing ratio 
distribution r (h) [Fig. 23], on the ray path geometry [Fig. 24], and 
on the refractivity distribution N(h) [Eq, 5]. Cumulative transmissivity 

along a ray path is defined by 

00 00 

A=Ja(A)dA and AT = jA<p{l)dZ, (43) 

h. h. 

1 1 

where djJ is the increment of the ray path and h. is the initial altitude. 

The atmosphere has only a weak dependence on geographic coordinates and 
in the first approximation it can be considered as a layered homogeneous 
medium which does not change with time. The numerical method used 
for the evaluation of equation (43) consisted of assuming a spherically 
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stratified atmosphere which was divided into slabs of thickness Ah. The 

altitude increments Ah were chosen so that each layer was quasi- 

homogeneous (Ap < 10 torr and/or AT <5“K across each slab) amounting 

to a number of n = 151 for h = 0 to 100 km (see Table 9). Attenuation a 

and dispersion A0 were calculated for each slab as described iu Section 

3. 1. 1. The clean, dry atmosphere sustains a constant mixing ratio up 

—6 

to h > 90 km for all gases which might have a measurable (r^ > 10 ) 

influence upon transfer properties (Fig. 23). An exception is ozone 
(and possible air pollutants), which is discussed briefly in Section 3. 2. 3 
The correct evaluation of a path trajectory would include refractive ray 
bending commonly expressed by Snell' s law (e. g. , Blake, 1972). The 
effect is most severe for a horizontal starting angle (6 = 90°) defining a 
tangential path and is zero for a vertical starting angle (9 = 0) defining 
a zenith path. We assume a straight line slant path with the increment 
di(9)« Aji(9) determined by the algorithm given in figure 24. Finally, 
a numerical integration based on Simpson' s rule was employed to com- 


pute the integrals of (43); i.e. , 

A 2 0. 5 (a +a .,) AX (9) 

n n n+1 n 

AT «=> 2 0. 5 (A0 + A0 ^,) AX (9) 

n n n+l n 


(44) 


A word of caution regarding these calculations should be noted 
at this point. Results obtained by the method described above can only 
qualify as analytical estimates allowing an overview on altitude - 
dependent transfer properties of the O^-MS, The following key -words 
underline the approximate nature: (a) empirical linewidth model, (b) 
straight line ray path, (c) numerical integration, (d) neglect of water 
vapor, (e) neglect of temporal variations of p and T, and (f) a more or 
less representative model atmosphere. The climate dependence of the 
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O^-MS is not too strong. The attenuation peaks of A , for example, 
decrease by about 10 percent in a tropical climate and increase by 
about ZO percent in an arctic climate when compared with results 
(Figs. 25, 17) for the U. S. Std. Atm. 62 (Reber et al. , 1970). Values 
for A and AT are, however, consistent in them selve s, and the spectra • 
could be scaled to results of reliable experiments. 

3. 2. 2 Zenith and Tangential Path Transmissivities Between 

49 and 72 GHz at Various Altitudes 

A computer program was developed for evaluations of equation 
(44). The routine was tested for the two extreme (zenith and tangential) 
cases of a slant path. We present examples of resvilts in the form of 
computer plots in Figures 25 through 39. The print -outs are too lengthy 
for reproduction here, but Table 11 gives an example of the format. 

Such listing can be made available upon request for any given set of 
variables. The path transmissions between a given initial altitude and 
outer space are obtained by stunming the contributions of 44 O^-MS 
lines, integrating over height, and finally integrating over frequency. 

For example, the computation of the results shown in figure 25 require 
for frequency increments of 20 MHz about 23 • 10^ individual calculations 
to cover the 49 to 72 GHz band which takes about 5 minutes on a CDC 
3800 machine. 

The zenith path spectra are labeled by the initial altitude h^ 
which was varied between h = 0 and 30 km in increments of 5 km. The 
tangential spectra are calcixlated for a straight-line path through the 
total atmosphere (satellite -to -sate Hite) and they are labeled by the 
height h^ , marking the closest approach of the tangent to sea level. 
Factors that determine tropospheric path loss are not only molecular 
attenuation but also stratification and path geometry which ceui introduce 
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multipath, diffraction, and ground reflection effects. To avoid con- 
sidering these problems, we limited ^ 5 km. The peak values of a 
zenith and a tangential path for h^ = h^ differ by about 1 to 30. The 
tangential spectra were extended to h^ = 50 and 7 5 km (Figs. 38, 39) to 
gain insight into the problem how deep a satellite -satellite link can 
penetrate the atmosphere. At these altitudes there is still appreciable 
attenuation at the line centers due to the long path length, but in between 
are many '’clear" channels with bandwidths on the order of 400 MHz. 

A few experimental results are available on zenith attenuation 
from h = 0 (Table 10). Carter et al. (1968) also reported some aircraft 
measurements for different altitudes between h^ = 2. 6 and 13. 7 km. 
Their data points for h^ = 9. 1 km are shown in figure 29. Figure 25 
illustrates the degree of agreement obtained between experiment and 
the CMR-Model analysis (Eqs. 32, 33) which is on the order of + 10 
percent. Table 10 gives the numerical values. 

Reber (1972) published recently more extensive data for Aon both 
sides of the O^-MS (Table 10). Special care was taken to eliminate 
water vapor attenuation, being typically on the order of < 1 dB total at 
zenith. The dry air data were fitted by yet another linewidth model 
(R-Model) which consists of a two-piece linear pressure dependence for 
the width y of the 60 GHz line complex excluding the 1 line -- 

Range 1, y < 52. 7 MHz, 


Range 2, y s 


y = 1. 88 p (300/ T) [MHz] 

52. 7 MHz (p^ = 20. 7 torr at T = 222 “K), 
y= 52. 7 + (1. 88/ 3) (p-p^) (300/ T) [MHz] . 


(45) 


We substituted the R-Model (Eq. 45) into the computer routine and, as 
expected, got a good fit to the Reber data. Both empirical width models 
(Eqs. 32, 45) are a somewhat arbitrary discription of an intensity dis- 
tribution which deviates from a linear addition of individual shapes. 
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The experimental zenith data (Table 10) strongly suggest overlap effects 
in the continuum spectrum {Sect, 2,4. 1), e specially when compared 
with A -values which have been calculated using the pre s sure -linear 
linewidth model (Eq. 19)* The experimental values, however, are only 
wing data and as such not really representative of the continuum 
spectrum. 

Center data cannot be obtained by the radiometric technique 
since the band between 55 and 65 GHz is practically opaque for any 
source in space. Reliable experimental information on the pressure 
dependence of the intensity distribution of the continuum spectrum 
can be obtained by a controlled laboratory experiment. An analysis of 
continuum data is based on the accurate knowledge of the line spec- 
trum and overlap effects as shown by the theoretical treatments of 
Mingelgrin (1972) and Dillon (1969)» 

The O^-MS offers between 55 and 65 GHz one feature not found 
at any lower frequency; namely, high attenuation and phase dispersion 
values related to the fairly stable dry part of the lower atmosphere. 

This affords shielding of satellite -satellite links against ground inter- 
ferences. The first 100 km of the atmosphere present a fairly com- 
plicated '^natural’* filter which can be used to advantage in special 
system applications. The strong frequency dependences of A (a) and 
AT(A0) permit frequency diversity techniques and signal preprocessing 
such as phase -keying. Figures 25 through 39 demonstrate that the 
atmosphere, which is opaque close to sea level, becomes transparent 
at higher altitudes in 16 channels each possessing a bandwidth of several 
hundred megahertz. One such channel is centered around 58, 8 GHz, and 
for a tangential path geometry the loss due to the O^-MS varies as 
follows 
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The curves presented in figures Z5 through 39 for cumulative trans- 
mission and zenith ray path properties are for one-way radio communi- 
cation cases. The distinction between the highly-structured line spec- 
trum (1) and the smooth continuum spectrum (II) can clearly be made 
within each cumulative spectrum. 

Aircraft can communicate with one another or with satellites 
between the O^-MS lines when flying above 10 kilometers while being 
shielded against the ground. Several atmospheric channels might be 
combined to increase the bandwidth by means of frequency transposition. 
The cumulative phase dispersion, AT , in the pass band centered at 58, 8 
GHz varies for a zenith path from -5. 7 rad per 350 MHz at h^ = 15 km, 
to -4. 5 rad per 680 MHz at h^= 30 km, 

3. 2. 3. Some Spectroscopic Properties of Ozone and Other 
Minor Atmospheric Gases. 

The mixing ratio for the main constitutents of clean, dry air is 
> -5 

constant up to h 90 km for all r ^ 10 (Fig, 23). Several minor 

K. 

constituents and gaseous pollutants exhibit, in addition to O and H O, 
spectral lines over the 40 to 140 GHz band (Table 12). These lines are 
part of rotational spectra due to molecular electric dipole moments. 
Linear rotor molecules (N^O, CO) exhibit an almost harmonic spacing 
of lines, while asymmetric rotors with three different principle moments 
of inertia (H^O, SO^, NO^, CH^O) display an irregular line pattern. 
Although the absorption coefficients of electric dipole transitions are 
considerably larger than those for these spectra are of no impor- 

tance to transfer properties of the lower (h< 20 km) atmosphere. 
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Trace concerLtratiorxs and strong foreign-gas -broadening essentially 
^Vash-out^^ these resonances. 

Ozone is practically the only molecule which causes some con- 
cern. The line shape is altered according to the manner in which the 
gas is distributed throughout the atmosphere. The increase of the 
mixing ratio with height produces line profiles of several megahertz 
width. Waters (1970) calculated 48 O^-lines between 39. 4Z6 and 
14Z. 172 GHz of which 14 lines are masked by the O^-MS. The 7 
strongest lines are listed in Table 12, and the maximum values for 
the total zenith attenuation indicate that the vicinity of these lines pro- 
bably should be avoided by ground -to -satellite links. On the other hand, 
one could tune -in' on these resonances and use phase or amplitude vari- 
ations of a transmitted signal as tracer for high-altitude ozone activity. 

3.2.4 Atmospheric Noise Due to Molecular Attenuation. 

Atmospheric noise reaching an antenna as a function of orienta- 
tion and frequency is, as derived from radiative transfer theory, an 
integral over the temperature T(X) and the attenuation rate a(j^) along a 
ray path (see Sect, 3. 2. 1). A ground-based antenna '*sees’^ the noise 
temperature (e . g. , Blake, 1972) 

00 J0 

= J T(i) a{i) exp. [-J a(s) ds] di [ °K] . {46) 

o o 

The tropospheric O^-MS attenuation between about 55 and 65 GHz is so 
high (opaque atmosphere), that most of the noise originates from the 
lowest layers. The noise temperature equals approximately the surface 
temperature since the integrand of (46) is large at i = 0 and rapidly ap- 
proaches zero for At the band edges of the -MS, the atmosphere be- 

comes semi-transparent, and T starts to drop when 30 dB (Eq. 43). 
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Altitude 



Relative Concentration by Volume . r 


Figure 23, Model profiles of atmospheric constituents showing their 
relative concentrations versus altitude up to h= 90 km 
(after Waters, 1970), denotes molecules with resonance 
spectra below 140 GHz (Tables 7, 12 and Sect, 2. 5). 
Numerical values see^ Table 1, 
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Figure 24, Algorithm used for the calculation of the differential slant 
path length A1 between two shells of a spherically strati- 
fied atmosphere. 6 - starting angle against zenith (n= 151 
for h=0 to 100 km). 
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Table 10. Zenith Oxygen Attenuation A from Sea Level 


F requency 
± 1 MHz 

a) 

Measu 

No. 

in 

Fig . 25 

red Values 

No. of 
Samples 



RMS- 
Error 
± 6A 

A 

Calculated Values 

R-^^ CMR-^^ Linear 

Model Model Model 

AAA 

[GHz] 





[%] 


[dB] 


50. 000 

W 

1957 

2 


15 

1.3 





R 

1972 


228 

7 . 1 

1 . 1*2 

1. 26 

1.55 

3 . 17 

50. 200 

R 

1 } 


130 

9.9 

1.31 

1.33 

1.64 

3 . 36 

50. 700 

R 

r 1 


120 

6.3 

1.58 

1 . 57 

1 . 92 

3 . 90 

51, 250 

R 

1 1 


226 

12 . 6 

1.98 

1.93 

2.35 

4.67 

51,750 

R 

r 1 

7 

638 

11 . 1 

2.43 

2.41 

2.91 

5 . 62 

52.220 

R 

r 1 


171 

4.8 

2.92 

3 . 07 

3.67 

6.81 

52.500 

SW 

1969 

6 

33 

^ZO 

5 . 1 

3.73 

4.37 

7 . 76 

53.400 

CMR 

1968 

1 

125 

5.2 

7 . 16 

6 . 88 

7 . 10 

12 . 18 


R 

1972 


123 

3.0 

7 . 06 




53.500 

SW 

1969 

6 

8 

«20 

9.3 





CMR 

1969 


no 

6.7 

7.82 

7 . 74 

7.82 

13 , 0 


R 

1972 


60 

4 . 0 

8.05 





R 

1972 


43 

2.9 

7.84 




53,811 

CMR 

1968 

1 

129 

5.4 

10.20 

9.37 

10 . 52 

15 . 33 


R 

1972 


120 

5.9 

9.59 




54.352 

CMR 

1968 

1 

26 

1,4 

15.39 : 





CMR 

1 r 

1 

57 

2 . 8 

15.62 ' 

14 . 53 

15.98 

21 . 17 


R 

1972 


44 

3.9 

13.91 




54.895 

CMR 

1968 

1 

50 

5 . 6 

23.35 

22 . 50 

24 . 35 

29 . 57 


R 

1972 


20 

4 . 1 

20.89 




66. 056 

R 

r 1 


122 

7.3 

8.61 

8 . 92 

10 . 25 

16 . 00 

66.565 

R 

11 


43 

0.5 

6.23 

6,28 

7.37 

12.46 

67, 016 

R 

11 


230 

6 . 3 

4.64 

4.79 

5 . 71 

10.21 

67.696 

R 

11 


144 

4.2 

3.08 

3.36 

4 . 09 

7. 82 

68. 144 

R 

11 

7 

402 

5 . 1 

2.58 

2 . 79 

3.42 

6 . 71 

68.680 

R 

1 1 


81 

3 . 0 

2.04 

2 . 31 

2 . 85 

5 . 70 

65 to 69 

TS 

1963 

3 

- 


18 to 2.8 


- 

- 

69.216 

R 

1972 


154 

13 . 5 

1.71 

1.97 

2.43 

4 . 92 

69.732 

R 

M 


81 

7.7 

1.46 

1.71 

2 . 12 

4 , 33 

70.0 

C 

1958 

4 

- 

16 

1.9 

1.62 

2 . 00 

4 . 07 

70.260 

R 

1972 

7 

447 

10 . 2 

1.37 

1.51 

1 . 88 

3.85 

72.0 

H 

1964 

5 

- 

25 

1.45 

1.08 

1 . 34 

2 . 75 


w 

“ Whitehurst et al, , 1957 

TS 

- Tolbert and Straiton, 1963 

R 

- Reber, 1972 

C 

- Coates, 1958 

SW 

CMR 

- Snider and Westwater, 1969 

- Carter et al. , 1968 

H 

- Hayes, 1964 


A survey of measured values of clear sky zenith attenuation for 
the frequency range 15 to 140 GHz has been made by Thompson III {1971; 
Fig. 9-11). 
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Captions for Figures Z5 to 

39 

— Computer Plots — 
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Table 12. Stronger Spectral Lines of Minor Atmospheric Gases 
in the 40 to 140 GHz Frequency Band 


Molecule C 

(type of 
rotor) 

lenter Freq. 

Max. Absorpt, 
Pure Gas, 300 ®K 

Zenith Absorpt. 
(max. ) 

A. 

Ground Level I 

Concentration 

clean/polluted 

Clec, Dipole 
Moment 

(asymmetr. ) 

(GHz J 

96.2288 

101.7368 

110.835 

124.086 

125.389 

136.883 
142 . 172 

LdB/kmJ 

426 ® ^ 
803 
1230 
1750 
1030 
300 
2330 

LdBj 

^ *) 
.16 > 

.3 

.4 

.5 

.75 

Lppm/ 

0 to 1 
0.07 j 

1 

3 to 1 

6 1 

at 1 

h= 1 
20 km ' 

vol J 1 

up to 
3 

[Debye or 
3 . 336-10 
Asm] 

0.53 


*T"Calculated by Waters ( 1970 ) for U. S. standard Atmosphere 
daytime ozone distribution (max. 5 x 10 molecules/cm a 

(1962) and 
t h = 20 ^) . 

NgO 

(linear) 

50.2460 

75.3696 

100.4917 

125.6137 

150.735 

20 ^ ^ 
65 ^^ 
152 
282 
48 o 


1 

0.5 1 

up to 
10 

0.167 

CO 

( linear) 

115.2712 

230.5380 

104 

838 


0.01 ’ 
to 

0.2_. 

up to 

150 

0.112 

SOg 

(asymmetr. 

12 lines listed between s 
53.529 and 131.015 CHz ' 


0 to 
0.02 

up to 
20 

1.615 

NOg 

( asymmetr. 

10 lines listed between 
) 40.358 and 41.278 GHz 


0 to 
0.01 

up to 
5 

0.316 

CHgO 

(asymmetr. 

8 lines listed between 
) 45.063 and 140.839 GHz 


0 to 
0.01 

up to 
[ 5 

2.34 

0 ^ (linear 

fine struG 
ture) 

^ (Table 7 

r 

) 1.56 

6.06 

120 (Fig. 25) 
»110 

209460 

Magnetic 

Moment 


a) Waters (1970) 

b) French (1968) 

c) Cord et al* (19b8) 
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Equation (46) was evaluated numerically for O and H O attenuation, for 
elevation angles 0 = 0® to 90® and h = 0 over the frequency band 0. 1 to 
100 GHz, by Blake (197Z, Fig. lOZ). The sensitivity of millimeter wave 
3^eceivers is not limited by atmospheric but by detection noise. Uncooled 
Schottky-diode mixers, for example, exhibit at 60 GHz presently single 
sideband noise temperatures on the order of 1Z00®K, 

In radio communication, data are transmitted through the atmos- 
phere in Such a way that the system channel parameters are matched to 
the propagation characteristics of the medium to maximize received 
data; in passive remote sensing, the medium itself provides the data, 
Meeks and Lilley (1963) introduced the idea of using O^-MS emission to 
remotely sense the temperature structure of the lower (h>30 km) 
atmosphere from an orbiting satellite. Their work has been extended 
by Lenoir (1968) and Groom (1971) and was applied by Westwater (1970) 
and Waters (1970, 1973) to ground-based radiometric sensing. All 
these treatments take advantage of microwave spectroscopic properties 
of and other atmospheric gases and their dependences upon frequency 
and meteorological conditions, which has been the subject of this report, 

F inding the temperature T(A) in equation (46) as a function of an 
independent variable [e.g, , p(h)] from observations of T at selected 
frequencies v requires numerical inversion methods. The reliability 
of the inversion depends to a great extent on the accuracy to which 
a(v, p, T, H, r^) due to the O^-MS is known, A general descrip- 
tion of the atmospheric O^-MS can be accomplished only by theory. 
However, existing data are insufficient for a proper test of the available 
theory. More laboratory work on the intensity distribution of the atmos- 
pheric O^-MS is required under careful control of conditions (Meeks a; 
Lilley, 1963; Wilheit, 1969: Liebe, 1969b; Waters, 1970). 
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4. CONCLUSIONS 

The review o£ molecular theory for millimeter wave properties 
of air, its translation into engineering terms (Section Z), the develop- 
ment of analytical schemes to predict propagation characteristics on the 
basis of mete or ological phenome na, and the description of computer 
routines together with many examples of unique atmospheric transfer 
properties (Section 3) shall provide a sound basis for the development 
of future systems operating in the 40 to 140 GHz frequency range. 

The main emphasis has been placed on the atmospheric O^-MS 
which dominates transfer properties throughout the 40 to 140 GHz band 
and existing gaps in the knowledge of the spectrum have been identified. 

It was attempted not to seriously violate real atmospheric situations; 
however, over the full range of meteorological variables, it remains to 
be established how reliable the discussed analytical schemes really 
are. Presently, the correlation between atmospheric attenuation and 
phase dispersion rates on one hand and meteorological conditions (Eq. 3) 
on the other hand is not very accurate. The correspondence between 
both sets of data needs to be improved by inputs from controlled labora- 
tory measurements on isolated line parameters and continuum spectra 
data. A manageable theory should be brought forward which predicts 
the continuum spectrum based on properties of the line spectrum, cir- 
cumventing more or less arbitrary empirical models. 
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